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I \- ABSTRACT 

I We derive the delay-time distribution (DTD) of type-la supernovae (SNe la) using a 

■ sample of 132 SNe la, discovered by the Sloan Digital Sky Survey II (SDSS2) among 
• ' 66,000 galaxies with spectral-based star-formation histories (SFHs). To recover the 

, best-fit DTD, the SFH of every individual galaxy is compared, using Poisson statis- 

'jp' tics, to the number of SNe that it hosted (zero or one), based on the method introduced 

Q I in Maoz et al. (2011). This SN sample differs from the SDSS2 SN la sample analyzed 

|H . by Brandt et al. (2010), using a related, but different, DTD recovery method. Fur- 

C/5 ' thermore, we use a simulation-based SN detection-efficiency function, and we apply a 

number of important corrections to the galaxy SFHs and SN la visibility times. The 
DTD that we find has 4ct detections in all three of its time bins: prompt (r < 0.42 Gyr), 

' intermediate (0.42 < r < 2.4 Gyr), and delayed (t > 2.4 Gyr), indicating a continuous 

^ , DTD, and it is among the most accurate and precise among recent DTD rcconstruc- 

l/^ ' tions. The best-fit power-law form to the recovered DTD is i-i-i2±o.08^ consistent with 

' generic ^ predictions of SN la progenitor models based on the gravitational- wave- 

, induced mergers of binary white dwarfs. The time integrated number of SNe la per 

O ■ formed stellar mass is Nsn/M = 0.00130 ± 0.00015 Mq-\ or about 4% of the stars 

\^ \ formed with initial masses in the 3 — 8 M0 range. This is lower than, but largely con- 

. sistcnt with, several recent DTD estimates based on SN rates in galaxy clusters and 

04 ' in local- volume galaxies, and is higher than, but consistent with Ns-m/M estimated by 

^ comparing volumetric SN la rates to cosmic SFH. 

^ , Key words: supernovae: general - methods: data analysis ~ galaxies: star formation 

'X: 

H ■ 



In the single degenerate (SD) model (Whelan & Iben 1974; 
;plodes in Type-la supernovae (SNe ,t ^ mooN u j r /tift^n 

, T . / Nomoto 1982), a carbon-oxygen white dwarr (WD) grows m 

tnp mV In nrno'pnitnr r^rnnlpm Itjpp 



elude the risk of systematic errors in cosmologicalinferences , , .., . iir-r^ , ■ , ,1 

° tational waves, with the more massive WD tidally disrupt- 

using SNe la as distance indicators, an incomplete picture of . , j_- 1 • ^ 1 ^ 

° , , r mg and accreting the lower-mass object, and at some point 

cosmic history (we do not know the physical identity 01 one 

Igniting. 

of the main metal producers), and ambiguity regarding the 

final outcome of two important stages in stellar and binary To date, neither the SD nor the DD model is clearly 

evolution - cataclysmic binaries and WD mergers. favored observationally. Both models, SD and DD, suffer 

from problems, theoretical and observational. One example 
is the lack of early-time UV shock signatures (Kasen 2010) 
* E-mail: maoz@astro.tau.ac.il due to ejecta hitting a SD red-giant donor (e.g., Nugent et al. 
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2011; Brown et al. 2012; Bloom et al. 2012), or lack of early 
radio and X-ray emission, expected if the SN blast wave 
were to encounter a circumstellar wind from a SD donor 
(e.g. Chomiuk et al. 2012, Horesh et al. 2012). Perhaps most 
dramatically, Schaefer & Pagnotta (2012) have shown that 
there is no remaining companion, down to luminosity limits 
corresponding to main-sequence K stars, in the remnant of 
a SN la from circa 1600 (confirmed as such with a light-echo 
spectrum by Rest et al. 2008). These observations disfavor 
non-degenerate donors. 

On the other hand, calcium, sodium, and other absorp- 
tion lines possibly associated with a SD donor star have been 
seen in some SN la spectra (Patat et al. 2007; Simon et al. 
2009; Sternberg et al. 2012; Dilday et al. 2012). Against the 
DD model, it has also long been argued that the merger of 
two unequal-mass WDs will lead to an accretion-induced col- 
lapse to a neutron star, i.e. a core-collapse SN (although not 
a typical one, due to the absence of an envelope), rather than 
a SN la (Nomoto & Iben 1985; Shen et al. 2012), but ways 
have been proposed to avoid this (e.g. Pakmor et al. 2010; 
Van Kerkwijk et al. 2010). Several "SN on hold" scenarios 
have also been proposed (Di Stefano et al. 2011; Justham 
2011; Kashi & Soker 2011; Ilkov & Soker 2012), where igni- 
tion is potentially long delayed. During this time the traces 
of the messy accretion process (or even of the donor itself) 
could disappear. More than one SN la channel could be at 
work, and this possibility has often been raised (e.g. Dilday 
et al. 2012). 

It has long been realized that SN la rates may discrimi- 
nate among progenitor models, via the "delay time distribu- 
tion" (DTD). The DTD is the hypothetical SN rate versus 
time that would follow a brief burst of star formation. It is 
the "impulse response" that embodies the physical informa- 
tion of the system, free of nuisances such as, in the present 
context, the star- formation histories (SFHs) of the galaxies 
hosting the SNe. It is directly linked to the stellar and bi- 
nary evolution timescales up to the explosion, so different 
progenitor scenarios predict different DTDs. Furthermore, 
apart from the progenitor problem, the SN DTD, together 
with SFH, dictates the heavy-element chemical evolution of 
individual galaxies and of the Universe as a whole. A well- 
determined DTD based on observations would be a valuable 
input to metal-enrichment history. Various theoretical DTDs 
have been proposed, whether based on detailed binary pop- 
ulation synthesis (e.g., Yungelson & Livio 2000; Nelemans 
et al. 2001; Han & Podsiadlowski 2004; Mennekens et al. 
2010; Ruiter et al. 2011; Meng et al. 2011; see Wang & Han 
2012, for a review), or on some simple physical considera- 
tions. In the DD model, the event rate ultimately depends 
on the loss of energy to gravitational radiation. If, at the 
end of the last common- envelope phase, the WD separation 
distribution follows a power law, then for a fairly large range 
of power-law indices around —1, the DTD for this model will 
also be a power-law in time, out to a Hubble time, ~ , with 
P not far from —1. Of course, the post-common-envelope 
separation distribution could be radically different, and not 
even be well-approximated by a power law, and thus a ~ 
DTD form should not be considered unavoidable (see, e.g., 
Ruiter et al. 2011). Nevertheless, the WD separation distri- 
bution that emerges from simulations does seem close to a 
power law of index —1 (see discussion in Maoz et al. 2012). 
In the SD models, in contrast, the DTD often cuts off after 



a few Gyr, due to a lack of evolved donors that can transfer 
enough mass to the WD. Here again, this expectation could 
be modified in the "on hold" scenarios mentioned above, 
which could introduce additional delays before explosion. 

Mannucci et al. (2005, 2006) and Sullivan et al. (2006), 
analyzing SN la rates as a function of galaxy colors, made 
the first direct estimates of SN la delay times, finding evi- 
dence for a range of delay times, with both "prompt" SNe 
la that explode within hundreds of Myr of star formation, 
and "delayed" SNe la that explode in populations of ages 
of at least a few Gyr. Totani et al. (2008) measured SN 
la rates in elliptical galaxies as a function of luminosity- 
weighted galaxy age, modeled by assuming a coeval stellar 
population, and deduced a DTD consistent with a form. 
Recently, a number of novel DTD recovery techniques, us- 
ing a variety of SN samples, environments, and redshifts, 
have been yielding mostly consistent DTDs. One approach 
to DTD recovery has been to measure SN la rates versus 
redshift in galaxy clusters, in which the observed SN la rate 
vs. cosmic time since the stellar formation epoch provides an 
almost direct measurement of the DTD. Another approach 
has been to measure the volumetric SN rate vs. redshift from 
field surveys. This rate will be the convolution of the DTD 
with the cosmic SFH. The cluster-based SN la DTD, as an- 
alyzed in Maoz et al. (2010), as well as the one emerging 
from volumetric SN la rates (Graur et al. 2011), also ap- 
pear to follow the t"^ form generically expected in the 
DD picture, although some controversies remain (see Maoz 
& Mannucci 2012). Very recent new and precise volumetric 
rates out to redshift z = 1 by Perrett et al. (2012) merge 
smoothly with the Graur et al. (2011) rates at higher z, and 
their analysis confirms the above DTD result. 

The above approaches to DTD recovery involve loss of 
information, by associating all of the SNe found at a given 
redshift with all the galaxies at that redshift. Brandt et al. 
(2010; BIO) and Maoz et al. (2011; Mil) simuffaneously 
introduced two closely related new methods for DTD re- 
construction that avoid this averaging. Rather than using 
SFHs averaged over many galaxies in some redshift bin, 
and comparing to the number of SNe found by a survey 
at those redshifts, the methods take into account the de- 
tailed SFH of each individual galaxy monitored by a SN 
survey. In the method used by Mil, the expected SN rate, 
for some assumed DTD, is compared to the observed SN 
rate in each galaxy. The DTD parameter space is searched 
for the DTD that best reproduces the observed individual 
galaxy SN rates. Mil applied the method to a subsample of 
the galaxies in the Lick Observatory SN Search (LOSS, Lea- 
man et al. 2011, Li et al. 2011a,b) that have SFHs based on 
Sloan Digital Sky Survey (SDSS; York et al. 2000) spectra, 
as recovered by Tojeiro et al. (2009), and to the SNe found by 
LOSS in those galaxies. Using a DTD with three time bins, 
prompt {t < 0.42 Gyr), "intermediate" (0.42 < t < 2.4 Gyr), 
and delayed {t > 2.4 Gyr), Mil found significant detections 
of both prompt and delayed SNe la, and an overall DTD 
consistent with a form. 

The DTD recovery method of BIO is similar to that 
of Mil in several respects. However, in BIO, rather than 
(as in Mil) comparing the SN numbers observed in each 
galaxy to the predictions of DTD models, the mean spec- 
trum of the SN host galaxies is compared to mean spectra 
formed with mock host-galaxy samples that result from each 
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of the tested DTD models. BIO applied their method to a 
sample of 77,000 galaxies in the "Stripe 82" region of the 
SDSS, having spectra with Tojeiro et al. (2009) SFH re- 
constructions, and to the SNe la that they hosted in the 
SDSS-II SN survey (SDSS2; Frieman et al. 2008; Sako et 
al. 2008). Like Mil, BIO detected, among their three SN la 
DTD time bins, significant signals in the t < 420 Myr and 
t > 2.4 Gyr bins. Furthermore, BIO separated their SNe la, 
based on light-curve shape, into "low-stretch" and "high- 
stretch" subsamples, and found that high-stretch SNe have 
a DTD signal mainly in the prompt bin, while low-stretch 
SNe have a DTD dominated by the delayed bin. 

In Maoz & Badenes (2010), the Mil method was ap- 
plied to a sample of SN remnants in the Magellanic Clouds, 
treating them as an effective SN survey, one in which the 
events are visible for tens of kyr. The SN rates in individual 
regions of the Clouds were compared to the detailed SFH 
of each region, as reconstructed by Harris & Zaritsky (2004, 
2009) using isochrone fitting to the resolved stellar popula- 
tions. In this sample, the limited number of SN remnants 
that originated from SNe la, of order 10, limited the preci- 
sion of the DTD, and hence allowed the significant detection 
only of the prompt SN la component. 

In the present paper, we apply the Mil DTD recov- 
ery method to a sample of SNe discovered in SDSS2. There 
are several differences in our analysis compared to the BIO 
analysis of the SDSS2 SNe. The DTD recovery approach, 
as described above, is different. We show, however, that the 
Mil method, when applied to the same sample used by BIO, 
and using the same assumptions as BIO, does give a DTD 
similar to the one found by BIO, indicating the approaches 
are equivalent. We construct a new SDSS2 SN la sample 
with careful attention to selection criteria, SN classification, 
and SN-host matching, resulting in a sample that is com- 
posed of SN events that are largely distinct from those in 
the BIO sample. With 132 SNe la, our sample is somewhat 
larger than the BIO sample, with its 101 SNe, and it is more 
complete, which is important in the DTD context. Finally, in 
our analysis, we use an empirical detection efficiency func- 
tion from the SDSS2 project, and we apply several neces- 
sary corrections, not included in BIO, to the calculation of 
the formed stellar mass in the monitored galaxies, and to 
the SN visibility time of each galaxy. These corrections do 
not affect much the form of the DTD, but do affect its nor- 
malization, resulting in a more robust time-integrated DTD. 
Further details on each of these new aspects of the present 
work are given in the body of this paper. 



2 DATA 

2.1 Galaxy sample 

The SDSS (York et al. 2000), in its first stage, was a survey 
of ~ 10^ deg^ of the north Galactic cap, consisting of imag- 
ing in five photometric bands {u, g,r,i, z), and 3"-aperture 
fibre spectroscopy of ~ 10® targets, mostly galaxies, with 
r < 18 mag. Tojeiro et al. (2009) performed spectral synthe- 
sis modeling of all galaxy spectra in the SDSS using their 
VESPA code (Tojeiro et al. 2007). VESPA uses aU of the 
available absorption features, as well as the shape of the 
continuum, to deconvolve the observed spectra and obtain 



an estimate of the SFH. The SFH of each galaxy consists of 
the amounts of stellar mass formed in fixed bins of look-back 
time. In order to recover the maximum amount of reliable 
information, the number of time bins used is variable, and 
depends on the quality of the data on each galaxy. At the 
highest resolution, VESPA uses 16 age bins, logarithmically 
spaced between 0.002 Gyr and to = 13.7 Gyr, the age of the 
Universe. When data do not have sufficiently high signal- 
to-noise ratio for a fully resolved reconstruction, pairs of 
adjacent time bins are averaged. This process may be re- 
peated down to the last two remaining bins. In the end, 
VESPA recovers each galaxy's SFH using a different set of 
time bins. 

VESPA masses are calculated assuming a Kroupa 
(2007) initial- mass function (IMF); all of our results will 
therefore include this assumption implicitly. Bell et al. 
(2003) have shown that the Kroupa IMF gives a similar 
total stellar mass to that of a "diet Salpeter" IMF, ob- 
tained by multiplying by 0.7 the total mass of the original 
Salpeter (1955) IMF, to account for the reduced number of 
low-mass stars. Thus, our results will be comparable to other 
SN rate studies, such as Mannucci et al. (2005, 2006), and 
other recent DTD reconstructions (Maoz et al. 2010; Maoz 
& Badenes 2010; Mil; Graur et al. 2011) that have assumed 
the diet-Salpeter IMF. 

During 3 years from 2005 to 2007, for 90 days a year, 
the SDSS2 repeatedly imaged, with a 4-day cadence. Stripe 
82, a 120° X 2.5° equatorial region, in the RA range from 
-60° to +60° (Frieman et al. 2008; Sako et al. 2008). As de- 
scribed in BIO, in Data Release 7 (DR7) of the SDSS there 
are about 77,000 Stripe-82 non-active galaxies with spec- 
tra that have VESPA SFH reconstructions by Tojeiro et al. 
(2009). Excluding galaxies beyond redshift z > 0.4 (beyond 
the detection limit of SNe la, see below) , several tens of ob- 
jects with blueshifts of about 10 km s~^ (likely stars that 
were misidentified as galaxies), and galaxies with VESPA 
fits with > 10, leaves about 67,000 galaxies. Finally, 
we consider only galaxies in the right ascension (RA) range 
-51° < aj20oo < +57° which Dilday et al. (2010; DIO) de- 
fined for the purpose of their SN rate measurements, and for 
which they performed detection-efficiency simulations (see 
below). This leaves 66,400 galaxies. About 600 of these are 
actually likely pairs of duplicate SDSS spectra, taken at dif- 
ferent dates, of the same galaxy (same coordinates to within 
3" and same redshifts to within 300 km s~^). The resulting 
overcount of the number of monitored galaxies by about 300 
(i.e. by a factor of 0.5%), is of no consequence to our anal- 
ysis. We consider these ~ 66, 000 galaxies as the spectral 
galaxy sample that was monitored for SNe in SDSS2. 

As in BIO and Mil, for the majority of the SDSS galax- 
ies, it is practical to separate the SFHs into no more than 
four time bins; 0-70 Myr, 70-420 Myr, 420 Myr-2.4 Gyr, 
and > 2.4 Gyr, corresponding to the bins labeled 24, 25, 26, 
and 27 in Tojeiro et al. (2009). Furthermore (also as in BIO 
and Mil), it is clear that VESPA often does not separate re- 
liably the SFHs in the first two bins, so we combine the first 
two time bins into a single bin, of — 420 Myr. We use the 
VESPA SFH reconstructions based on the Maraston (2005) 
spectral synthesis models with a single dust component. 

We scale down the stellar masses listed in the VESPA 
database by a factor of 0.55, for the following reasons. To 
account for the limited 3" SDSS fibre aperture, Tojeiro et al. 
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(2009) scaled up the formed masses reconstructed from the 
SDSS spectrum of each galaxy, based on the difference be- 
tween the z-band magnitudes for the central 3" and the total 
Petrosian magnitudes (Blanton et al. 2001), where both of 
these magnitudes are based on the SDSS imaging. However, 
the spectral magnitudes, obtained by integrating the SDSS 
spectrum over the relevant filter curves, and which matched 
the fibre magnitudes in SDSS Data Release 5 (DR5), no 
longer do so in the later data releases, including DR7 on 
whose spectra the reconstructed VESPA masses are based 
(R. Tojeiro 2010, private communication). This recalibration 
leads (Adelman-McCarthy et al. 2008) to a 0.35 mag shift (a 
factor of 0.72). However, a further factor of 0.75 to 0.80 is re- 
quired (R. Tojeiro 2010, private communication) in order to 
make the VESPA galaxy masses match the mass estimates 
of the same galaxies based on pre-DR6 SDSS spectroscopic 
calibrations (Gallazzi et al. 2005; Brinchman et al. 2003). 
Apparently, this is the result of the DR7 spectral calibration 
being, on average, slightly redder than pre-DR6 spectra at 
wavelengths > 6500 A. This, in turn leads to higher mass 
estimates in DR7. Furthermore, we have found a similar fac- 
tor 0.75 scaling by comparing the VESPA masses available 
for 1282 galaxies in the sample of Mannucci et al. (2005) 
to the masses of those galaxies based on their B and K- 
band photometry, and using the relations of Bell & de Jong 
(2001). Finally, a similar factor of 0.75 is suggested by com- 
paring VESPA masses and PEGASE (Rocca-Volmerange et 
al. 2011) masses, obtained for a sample of 13 galaxies by 
Neill et al. (2009), by fitting their spectral energy distri- 
butions to broad-band photometry, ft is possible that, in 
fact, it is the DR7 spectra and their VESPA reconstructions 
that give the more accurate stellar masses, in which case it 
would be the previous SN rates per unit mass and DTDs 
(e.g. Mannucci et al. 2005; Li et al. 2011; Mil) that would 
have to be scaled down by a factor 0.75. As this is still un- 
clear, we have chosen here to apply this additional factor of 
0.75, resulting in a full factor 0.55 correction, to the Stripe 
82 galaxy VESPA SFHs. The full factor of 0.55 was also ap- 
plied to the galaxy SFHs in Mil, which were also based on 
DR7 spectra. Neither part of this correction, however, was 
applied in BIO, and as a result the DTD values in BIO will 
be about a factor of 2 lower than here and in Mil, just for 
this reason. 

The VESPA galaxies in Stripe 82 are our chosen sample 
of galaxies monitored for SNe. For each galaxy, we therefore 
need to know the SN la detection efficiency and the effective 
visibility time. Dilday et al. (2008, 2010) have carried out 
detailed efficiency simulations for the purpose of rate mea- 
surements, using fake SNe planted in the SDSS2 images in 
the course of the survey. These images were then searched 
for SN candidates, and subjected to detection and classifi- 
cation criteria, just like the real SNe. We adopt the fidu- 
cial SN la detection efficiency, as a function of redshift, as 
given in figure 8 of DIO. We approximate this function with 
a piecewise linear dependence - a constant for z < 0.175: 
e(z < 0.175) = 0.72; and a linear decrease to e = be- 
tween z = 0.175 and z = 0.4: e{z) = -3.2z + 1.28. This 
approximation matches the empirical curve to a few percent 
at all redshifts. The < 100% efficiency, even at low redshifts 
where all SNe la are detected by the survey, mainly reflects 
the selection criteria of the survey, which required that SNe 
included in the final sample be observed near maximum light 



and also well after maximum light. This effectively rejects 
SNe that exploded very near the beginning and near the end 
of every observing season (see Dilday et al. 2008, and fur- 
ther below). This detection efficiency function differs from 
the one used in BIO. In BIO, a particular functional form for 
e{z) was assumed, with a number of free parameters. These 
parameters were then tuned such that mock samples of SNe 
la exploding in the sample galaxies, for an assumed DTD, 
would have a similar redshift distribution to that of the real 
SN la sample. Overall, the BIO efficiency function assumes a 
lower assumed detection efficiency, which leads to increased 
DTD values. 

SDSS2 was a "rolling survey", in which the observer- 
frame visibility time for each monitored galaxy was, in prin- 
ciple, just the the duration of the survey, 269 days (DIO). 
The loss of visibility time at the "edges" of each observing 
season is already accounted for in the detection efficiency 
function, as discussed above. In the rest frame of a galaxy 
at redshift z, the SN visibility time is reduced by a factor 
1-1-2 due to cosmological time dilation. This reduction was 
inadvertently omitted in BIO (although the tuning of the 
detection efficiency parameters in order the match the SN 
redshift distribution may have partly accounted for this ef- 
fect). 

2.2 Supernova sample 

We turn next to the SN sample. About one-thousand likely 
SN candidates were discovered in the course of SDSS2, and 
some of these candidates were then followed up with spectro- 
scopic observations. We select our SN la sample from among 
two subsamples defined by DIO and one subsample defined 
by Sako et al. (2011; Sll). From these three subsamples, we 
associate between SNe and Stripe 82 galaxy hosts. DIO im- 
posed selection criteria for inclusion of SN la candidates in 
their samples, criteria that were applied also in their detec- 
tion and classification simulations and are therefore reflected 
in their detection efficiency function (see above). They in- 
cluded: 

1. At least one photometric measurement with signal-to- 
noise ratio S/N > 5 in each of three bands, g, r, and i; 

2. A photometric observation at least 2 days (in the SN rest 
frame) before maximum of the best-fit light-curve; 

3. A photometric observation at least 10 days (rest frame) 
after maximum of the best-fit light-curve; 

4. A good fit to a SN la light curve based on the MLCS2k2 
program (Jha et al. 2007), with fit criteria detailed in DIO; 
and 

5. RA in the range —51° < aj2ooo < +57°. 

With these criteria, DIO obtained a subsample of 312 spec- 
troscopically confirmed SNe la, listed in their table 2. A 
second subsample, consisting of an additional 148 likely SNe 
la satisfying the above criteria, was obtained by DIO (and 
listed in their table 3) by collecting events that do not have 
spectroscopic confirmation as SNe la, but are very likely SNe 
la based on their photometric light curve analysis (criterion 
4, above), and which are associated with a host galaxy at a 
spectroscopic redshift consistent with the photometric red- 
shift from the light curve analysis. 

Sll also assembled a "photometric" sample of 210 
SDSS2 SNe la (listed in their table 3), lacking spectroscopic 
confirmation of the SN but with host galaxies having spec- 
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troscopic redshifts. The SN la classifications of Sll were 
based on photometric criteria that are not identical, but 
quite similar to those of DIO: 

1. At least one photometric measurement with signal-to- 
noise ratio S/N > 5 in two of three bands, g, r, and i; 

2. A photometric observation at —5 < t < +5 days (rest 
frame) around maximum of the best-fit light-curve; 

3. A photometric observation at+5<t<-|-15 days (rest 
frame) after maximum of the best-fit light-curve; and 

4. A good fit to a SN la fight curve based on the PSNID 
program (Sako et al. 2008), with criteria detailed in Sll. 
Applying the DIO criterion on RA (No. 5, above) to the Sll 
sample reduces it from 210 to 204 SNe. There are 107 SNe 
common to these DIO and Sll photometric SN la samples. 

We have searched for associations between the Stripe 
82 VESPA galaxies and the SNe la from these three sub- 
samples (the DIG spectroscopically confirmed sample, and 
the DIO and Sll photometric samples), by requiring a pro- 
jected physical SN-host separation of < 30 kpc, and a differ- 
ence of Az < 0.001 between the spectroscopic galaxy red- 
shift and the SN redshift (spectroscopic for the confirmed 
sample, photometric, based on the light-curve fitting, for 
the two other samples). Among the matches found in prac- 
tice, except for seven cases, the redshift difference is actually 
Az < 0.0005. The SN-host separation is always < 26 kpc, 
and except for 11 cases it is < 20 kpc. We find 61 galaxy 
matches to the spectrally confirmed SN la sample. All but 
one of these events, the SN with SDSS designation 15467, 
have lAU designations. We further match 32 SNe from the 
DIO photometric sample, and 62 SNe from the Sll photo- 
metric sample, after imposing on the latter the DIO RA cri- 
teria (which removes 3 SNe). Of the photometric SNe, 23 are 
common to the two samples, resulting in 71 unique photo- 
metric SNe. Together with the spectroscopic SN la sample, 
we therefore have a final sample of 132 SNe la that were 
hosted by the Stripe 82 VESPA galaxies (in the restricted 
RA range) during the SDSS2 SN survey, and were detected 
and classified with the efficiency estimated by DIO. Con- 
tamination of the photometric samples by non-SN la events 
is estimated at 3% (DIO) and 6% (Sll). We therefore ex- 
pect that only three or so of the 132 SNe are misclassified 
core-collapse SNe. Table 1 lists the adopted SN la sample. 

The SN la sample assembled by BIO in their DTD re- 
covery analysis is similar in size (101 SNe) to the present 
sample (132 SNe), but quite different in composition. BIO 
selected from the SN list available in the SDSS2 website 
the events that were listed as spectroscopically confirmed 
SNe la with lAU designations, and with light curves of suffi- 
cient quality for the light-curve analysis done by BIO. How- 
ever, this SN compilation is not complete. The DIO sam- 
ple of spectrally confirmed SNe has 12 SNe with Stripe-82 
host-galaxy matches that are not included in the BIO sam- 
ple. More importantly, the BIO sample omits the many SNe 
that were not observed spectroscopically, but are nonetheless 
bona- fide SNe la, as confirmed by their photometric analysis. 
On the other hand, BIO did include four SNe from the 2004 
observing season, before the full SDSS2 SN survey. In the 
end, only 50 SNe are common to BIO and to the present sam- 
ple. Finally, as already noted, this compilation of SNe did 
not have an associated detection and classification efficiency 
function, which BIO therefore estimated by reproducing the 
SN host redshift distribution. 



3 DTD RECOVERY METHOD 

We briefly repeat here the description of the DTD recovery 
method of Mil. For a sample of A'^ galaxies, the SN rate 
in the ith galaxy observed at cosmic time t is given by the 
convolution 

U{t) = r S,{t-T)^{r)dT, (1) 
Jo 

where: Si{t) is the star- formation rate versus cosmic time, of 
the ith galaxy (stellar mass formed per unit time); ^'(r) is 
the DTD (SNe per unit time per unit stellar mass formed); 
and the integration is from the Big Bang {t — 0) to the 
time of observation. We assume that the DTD is a universal 
function: it is the same in all galaxies, independent of en- 
vironment, metallicity, and cosmic time — neglecting such 
possible dependences (e.g., variations in the IMF with cos- 
mic times or with environment would also lead to a vari- 
able DTD). The DTD is recovered by inverting a linear, 
discretized version of Eq. [T] The SFHs of the i — 1, 2, TV 
galaxies monitored as part of a SN survey are known with 
a temporal resolution that permits binning the stellar mass 
formed in each galaxy into j = 1,2, ..,K discrete (but not 
necessarily equal) time bins, where increasing j corresponds 
to increasing lookback time. For the ith galaxy in the survey, 
the stellar mass formed in the jth time bin is mij . The mean 
of the DTD over the jth bin (corresponding to a delay range 
equal to the lookback-time range of the jth bin in the SFH) 
is ^'j . Then the integration in Eq. [1] can be approximated 
as a sum, 

K 

ri ^m,,*,, (2) 

where Vi, the SN rate in a given galaxy, is measured at a 
particular cosmic time, corresponding to the galaxy redshift, 
Zi. Given a survey of A'' galaxies, each with an observed 
SN rate, r;, and a known binned SFH, rriij, one could, in 
principle, algebraically invert this set of linear equations and 
recover the best-fit parameters describing the binned DTD: 
* = (*i,*2,...,*x). 

In practice, on human timescales SNe in a given galaxy 
are rare events (vi 1 yr"'^). Supernova surveys therefore 
monitor many galaxies. For a given model DTD, ^, the ith 
galaxy will have an expected number of SNe 

\i = riti^i, (3) 

where tiei is the effective visibility time during which a SN 
of a particular type would have been visible. In the present 
context of the SDSS2 rolling survey, ti — 269 d/{l + Zi), and 
is the detection efficiency based on the DIO simulations 
(which take into account the actual on-target monitoring 
time, the SN detection limits at the redshift of the galaxy, 
the light-curve selection criteria, and the detection and clas- 
sification efficiencies). Since Xi <^ 1, the number of SNe 
observed in the ith galaxy, Ui, obeys a Poisson probability 
distribution with expectation value Xi, 

P{n,\X,)^{e-^'X'^')/n,\, (4) 

where n; is for most of the galaxies, and 1 for some of the 
galaxies. 

To deal with this Poisson aspect of the problem. Mil 
introduced a non-parametric, maximum-likelihood method 
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Table 1. Supernova and host galaxy sample 



Supernova 


Host galaxy 




SN-host Sep. 


lAU 


ID RA Dec z 


RA Dec 


z 


" kpc 


ID 


(1) (2) (3) (4) 


(5) (6) 


(7) 


(8) (9) 


(10) 


Spectroscopically confirmed SNe 



762 


15, 


,53518 


-0, 


,87907 


0, 


1910 


15, 


,53606 


-0, 


,87966 


0, 


19154 


3, 


,82 


12, 


,18 


2005eg 


1032 


46, 


,79565 


1, 


,11952 


0, 


1300 


46, 


,79592 


1, 


,12000 


0, 


12978 


1, 


,97 


4, 


,56 


2005ez 


1112 


339, 


,01761 


-0, 


,37527 


0, 


,2580 


339, 


,01691 


-0, 


,37489 


0, 


,25777 


2, 


,87 


11, 


,47 


2005fg 


1371 


349, 


,37375 


0, 


,42929 


0, 


1190 


349, 


,37369 


0, 


,42966 


0, 


,11906 


1, 


,36 


2, 


,92 


2005fli 


2561 


46, 


,34335 


0, 


,85829 


0, 


1180 


46, 


,34430 


0, 


,85972 


0, 


,11841 


6, 


,19 


13, 


,23 


2005fv 


2689 


24, 


,90027 


-0, 


,75879 


0, 


1620 


24, 


,90004 


-0, 


,75796 


0, 


16153 


3, 


,10 


8, 


,62 


2005fa 


2992 


55, 


,49692 


-0, 


,78269 


0, 


1270 


55, 


,49731 


-0, 


,78294 


0, 


12656 


1, 


,67 


3, 


,79 


2005gp 


3241 


312, 


,65143 


-0, 


,35416 


0, 


,2590 


312, 


,65280 


-0, 


,35423 


0, 


,25900 


4, 


,95 


19, 


,87 


2005gh 


3592 


19, 


,05240 


0, 


,79183 


0, 


,0870 


19, 


,05289 


0, 


,79061 


0, 


,08661 


4, 


,74 


7, 


,70 


2005gb 


3901 


14, 


,85039 


0, 


,00252 


0, 


,0630 


14, 


,85049 


0, 


,00266 


0, 


,06286 


0, 


,61 


0, 


,74 


2005ho 


6057 


52, 


,55355 


-0, 


,97467 


0, 


,0670 


52, 


,55368 


-0, 


,97448 


0, 


,06706 


0, 


,83 


1, 


,06 


2005if 


6295 


23, 


,67295 


-0, 


,60544 


0, 


,0800 


23, 


,67429 


-0, 


,60420 


0, 


,07961 


6, 


,58 


9, 


,89 


2005js 


6406 


46, 


,08860 


-1, 


,06301 


0, 


1250 


46, 


,08863 


-1, 


,06312 


0, 


12463 


0, 


,41 


0, 


,91 


2005ij 


6558 


21, 


,70165 


-1, 


,23814 


0, 


,0570 


21, 


,70190 


-1, 


,23815 


0, 


,05740 


0, 


,90 


1, 


,00 


2005hj 


7876 


19, 


,18237 


0, 


,79452 


0, 


,0760 


19, 


,18279 


0, 


,79361 


0, 


,07637 


3, 


,59 


5, 


,20 


2005ir 


12780 


322, 


15454 


1, 


,22804 


0, 


,0490 


322, 


,15671 


1, 


,23017 


0, 


,04944 


10, 


,94 


10, 


,58 


2006eq 


12856 


332, 


,86542 


0, 


,75589 


0, 


,1720 


332, 


,86539 


0, 


,75559 


0, 


,17171 


1, 


,09 


3, 


,17 


200611 


12874 


353, 


,96448 


-0, 


,17724 


0, 


,2450 


353, 


,96329 


-0, 


,17659 


0, 


,24502 


4, 


,89 


18, 


,83 


2006fb 


12950 


351, 


,66742 


-0, 


,84032 


0, 


,0830 


351, 


,66730 


-0, 


,84061 


0, 


,08271 


1, 


,13 


1, 


,76 


2006fy 


13070 


357, 


,78500 


-0, 


,74645 


0, 


1990 


357, 


,78491 


-0, 


,74657 


0, 


19855 


0, 


,56 


1, 


,83 


2006fu 


13072 


334, 


,95929 


0, 


,02439 


0, 


,2310 


334, 


,96069 


0, 


,02368 


0, 


,23062 


5, 


,67 


20, 


,86 


200611 


13354 


27, 


,56474 


-0, 


,88735 


0, 


1580 


27, 


,56472 


-0, 


,88671 


0, 


,15759 


2, 


,30 


6, 


,26 


2006hr 


13511 


40, 


,61224 


-0, 


,79419 


0, 


,2380 


40, 


,61127 


-0, 


,79422 


0, 


,23758 


3, 


,49 


13, 


,14 


2006hh 


14279 


18, 


,48869 


0, 


,37156 


0, 


,0450 


18, 


,48826 


0, 


,37143 


0, 


,04543 


1, 


,62 


1, 


,45 


2006hx 


14284 


49, 


,04921 


-0, 


,60105 


0, 
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49, 


,04937 


-0, 


,60099 


0, 


18109 


0, 


,61 


1, 
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14377 


48, 
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-0, 
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0, 


1390 
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0, 
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1, 
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4, 


,43 
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-0, 
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0, 
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0, 
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0, 
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3, 
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0, 
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0, 
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1, 
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5, 
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2006ju 
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0, 
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0, 
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0, 
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0, 
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0, 
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2, 
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2, 
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0, 
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2, 
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0, 
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0, 
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3, 
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0, 
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1, 
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4, 
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4, 
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,71829 
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0, 
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,71881 
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,19581 


0, 


,17499 


3, 
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0, 
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0, 
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0, 
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0, 
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4, 
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1, 
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1, 
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2, 
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0, 
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1, 
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2, 
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8, 
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0, 
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0, 
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6, 
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9, 
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0, 
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0, 
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0, 
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3, 
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9, 
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0, 
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0, 
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0, 
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6, 
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0, 
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0, 
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0, 
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2, 
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11, 
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0, 
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1, 
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3, 


,78 


2007om 


20084 
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0, 


,1300 


312, 


.80679 
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,12946 


7, 


.74 


17, 


,87 


2007ph 
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SDSS2 Supernova Delay Time Distribution 



Supernova Host galaxy SN-host sep. lAU 
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Supernova Host galaxy SN-host sop. lAU 



ID 


RA 


Doc 


z 


RA 


Dec 


z 


II 


kpc 


ID 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


20331 


7.53944 


1.24604 


0.1845 


7.53964 


1.24591 


0.18449 


0.84 


2.59 




20480 


357.27420 


0.91829 


0.1678 


357.27429 


0.91820 


0.16776 


0.47 


1.35 




20626 


8.47575 


-0.59308 


0.2760 


8.47590 


-0.59297 


0.27628 


0.67 


2.80 




20721 


323.18481 


-0.62219 


0.2120 


323.18481 


-0.62275 


0.21183 


2.01 


6.94 




20726 


42.27992 


-0.15946 


0.3201 


42.27989 


-0.15947 


0.32015 


0.13 


0.62 




20787 


50.51086 


-0.44245 


0.2707 


50.51080 


-0.44254 


0.27066 


0.40 


1.67 




20788 


51.67339 


-0.47783 


0.3940 


51.67345 


-0.47777 


0.39347 


0.31 


1.67 




21709 


326.52615 


-1.04962 


0.1590 


326.52771 


-1.04932 


0.15900 


5.71 


15.66 




21872 


7.61328 


-0.59137 


0.2286 


7.61338 


-0.59148 


0.22854 


0.52 


1.91 




22006 


48.78285 


-0.96254 


0.3956 


48.78292 


-0.96260 


0.39556 


0.32 


1.72 





Column header explanations: 
(1)- SDSS-II supernova identifier; 

(2-4) - Supernova right ascension and declination, J2000, in decimal degrees, and redshift; 

(5-7) - Associated host galaxy right ascension and declination, J2000, in decimal degrees, and redshift; 

(8-9) - Supernova-host separation, in arcseconds, and in projected kpc (assuming cosmological parameters h = 0.7, 

nm = 0.3, Ha = 0.7); 

(10) - Supernova lAU designation (for spectroscopically confirmed SNc only). 



to recover the DTD and its uncertainties, based on Cash 
(1979) statistics. Considering a set of model DTDs, the like- 
hhood of a particular DTD, given the set of measurements 
ni, njv, is 

JV 

L = '[[P{n,\\i). (5) 

i = l 

More conveniently, the log of the likelihood is 

N N 

lnL = ^lnP(ni|AO = -^A, -^ ^ InA, (6) 

i — l i — 1 ^i — 1 

where only galaxies hosting SNe contribute to the second 
term. The best-fitting model can be found by scanning the 
parameter space of the vector ^' for the value that maxi- 
mizes the log-likelihood. This procedure naturally allows re- 
stricting the DTD to have only positive values, as physically 
required (a negative SN rate is meaningless). 

The covariance matrix Cjk of the uncertainties in the 
best-fit parameters can be found (e.g.. Press et al. 1992) by 
calculating the curvature matrix, 

^ 1 dHnL _ A d[lnP{ni\X,)] 9[lnP(w,|A.)] 



JV 

= ^t?(n,/A, - l)^mijmife, (7) 
1=1 

and inverting it, 

[C] = [a]-\ (8) 

We have applied this DTD recovery method to the sam- 
ple of monitored galaxies, and the SNe la that they hosted 
in SDSS2, as presented in Section (2] 



4 RESULTS 

4.1 Delay-time distribution 

Figure [1] shows the recovered DTD for our sample of 132 
SNe la discovered during the SDSS2 among the 66,000 
Stripe-82 galaxies that have VESPA SFHs. The numeri- 
cal values and uncertainties of the DTD are listed in Ta- 
ble 2. Several aspects of the DTD are notable. First, we 
detect, at > 4(7 significance, a DTD signal in all three 
VESPA time bins, including the intermediate time bin, 
0.42 < t < 2.4 Gyr. BIO and Mil found highly signifi- 
cant detections in the prompt bin and in the delayed bin, 
but only ~ 2a DTD detections at intermediate delays (al- 
beit consistent with the level expected from a t^^ DTD). 
Our new result provides the strongest support to date for 
previous indications for a continuous DTD (Mannucci et al. 
2006; Totani et al. 2008; Maoz et al. 2010). These previous 
DTD reconstructions, however, in place of the actual SFHs, 
used luminosity-weighted galaxy ages (Totani et al. 2008) 
or simple assumed SFHs (Mannucci et al. 2006, Maoz et al. 
2010). The present result is the most significant and precise 
detection of a DTD signal in all three of the time-delay bins 
defined here, in a DTD reconstruction that uses individual 
galaxy SFHs. 

Also shown in Fig. [TJ for comparison, are a number 
of other recent DTD estimates. The value we find for the 
prompt DTD component, 'i'l, is in excellent agreement with 
the one found in Mil. The values of the intermediate and de- 
layed components, 'I'2 and ^'3, respectively, are a factor ~ 2 
lower than found by Mil. In the case of ^'2, the difference is 
not statistically significant, and furthermore the Mil value 
may have been overestimated due to some residual "leak" of 
signal from ^'i (see Mil, and below). The difference in ^3 
values between Mil and the present work, however, is more 
significant, at the ~ 2a level. Furthermore, both the current 
and the Mil ^'3 measurements are lower than the DTD val- 
ues implied by a number of independent measurements of 
SN la rates in galaxy clusters, as compiled in Maoz et al. 
(2010), and also shown in Fig.[T] On the other hand, as also 
shown in the figure, our new ^3 measurement is consistent 
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Figure 1. Recovered SN la DTD for the SDSS2 SN sample. Filled red circles mark the best- fit DTD values for each time bin, whose 
time range is indicated by the horizontal error bars. Vertical error bars show the Gaussian la uncertainties. Red dashed line is the best 
power-law fit, with index —1.12, to the recovered DTD. Previous DTD measurements also shown are: Mil analysis of the nearby LOSS 
sample (empty blue squares, with central values slightly shifted to the left, for clarity); Maoz et al. (2010) analysis based on SN la rates 
in galaxy clusters (empty black circles); and power law DTD found by Graur ct al. (2011), based on comparison of volumetric SN 
la rates and cosmic SFH (solid green line, with dotted lines marking the Itr range). 



in level with the DTD level at long delays found by Graur 
et al. (2011) from comparison of volumetric SN la rates to 
cosmic SFH. Our result thus revives previous speculation 
that the SN la rate in old populations may be enhanced by 
a factor of a few in cluster environments, compared to field 
environments (Sharon et al. 2007; Mannucci et al. 2008). The 
intermediate ^3 value found in Mil could be due to a rela- 
tively large fraction of LOSS galaxies in clusters. Among all 
the LOSS galaxies, 20% are in galaxy clusters that are listed 
in the NED database, based on the same criteria as applied 
in Mannucci et al. (2008; projected separation < 1.5 Mpc, 
redshift difference < 1000 km s~^). The fraction of the early- 
type LOSS galaxies in clusters will be higher. Assuming the 
SDSS2 galaxy sample, in contrast, is dominated by a field 
population, a factor-few enhancement of the SN la rate in 
clusters could lead to the higher ^'3 value measured in Mil, 
compared to here. We defer to future work a more detailed 
examination of this possibility regarding the LOSS sample. 

It is important to note that the errors shown in Fig. [T] 
are statistical only. A systematic source of error, one that 
is difficult to estimate, is associated with the SFH recon- 
struction of each galaxy (see Tojeiro et al. 2009, Mil, and 
Section 4.3, below). 

We have used minimization to fit model power-law 
functions, '^{t) oc t*^, to the recovered binned DTD, with the 
amplitude and the power-law index as free parameters. The 
model is averaged over each time bin, to obtain the value 
that is compared to the observed DTD value for that bin. 
The best-fitting power-law has an index of /? = — 1.12±0.08. 



Our result augments a growing number of indications for a 
DTD that is a power law with a slope of ~ —1. However, 
Fig. [T] suggests that a better description may be a function 
that steepens somewhat at long delays, compared to a sin- 
gle power law. Further studying the DTD at such a level of 
detail will require better temporal resolution. The present 
three-bin time resolution is dictated by the spectral resolu- 
tion and the signal-to-noise ratio of the SDSS spectra. 

4.2 SN la production efficiency 

As discussed in Maoz & Mannuci (2012), and references 
therein, another interesting observable is the DTD normal- 
ization. It can be discussed in the context of the integral of 
the DTD over the age of the universe, to, 

Nsn/M = / " *(t) dt, (9) 
Jo 

which gives the total number of SNe that eventually explode, 
per unit stellar mass formed in a short burst of star forma- 
tion. The integral in Eq. [9] can be approximated by a sum 
over the binned DTD, 

Nsn/M^ J2 ^j^tj- (10) 

From the simulations described in Mil, for declining power- 
law DTDs represented with three time bins, Eg.llOltvpicallv 
underestimates Nsn/M systematically by ~ 10 — 20%. 

For the SDSS2 sample, we obtain Ni^/M = 0.00130 ± 
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Table 2. Recovered Delay Time Distribution 



Sample 




A^SN 


*i 


^2 


^,3 


Nsk/M 








0-0.42 Gyr 


0.42-2.4 Gyr 


2.4-14 Gyr 




(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


Pull 


66,000 


132 


140 ± 30 


25.1 ±6.3 


1.83 ± 0.42 


0.00130 ±0.00015 


High s 


66,000 


67 


87± 24 


10.5 ±4.2 


0.78 ±0.32 


0.00066 ± 0.00011 


Low s 


66,000 


45 


24 ±15 


5.7 ±3.6 


0.96 ± 0.27 


0.00032 ± 0.00007 



Column header explanations: 

(1) - Sample used to derive DTD: Full sample, high-stretch sample, or low-stretch sample. 

(2) - Number of galaxies in sample. 

(3) - Number of SNe in sample. 

(4-6) - DTD rates and 68% uncertainty ranges, in units of IQ-l'* SNe yr-^ M~^. 
(7) - Time-integrated DTD, in units of SNe M~^. 



0.00015 MgV This is lower than the Ni^/M = 0.0020 ± 
0.0006 Mq found in Mil for the LOSS sample, although 
consistent with it, given the larger uncertainty of the Mil 
result. The present result is also consistent with a lower es- 
timate of iVia/Af = 0.0010 ± 0.0005 Mq^ found by Graur 
et al. (2011), from comparison of volumetric SN la rate 
evolution and cosmic SFH. A very recent analysis of the 
volumetric rates by Perrett et al. (2012) gives Nia./M ^ 
0.0006 - 0.0010 Mg^ (after converting their result to a 
diet-Salpeter IMF), similar to the Graur et al. (2011) re- 
sult. As was the case when considering the DTD amplitude 
^3, above, the estimates of Nii^/M from galaxy clusters are 
higher when based on cluster SN la rates vs. redshift, and 
even more so when based on cluster iron mass content (Maoz 
et al. 2010; see Maoz & Mannucci 2012, table 1). For exam- 
ple, the 2(7 lower limit is Ni^/M > 0.0017 from SN la 
rates in clusters, and Nis./M > 0.0034 Mg^ from the iron 
mass content. It remains to be seen whether this is the re- 
sult of some systematic error, or a real enhancement of SN 
la production in cluster environments (Sharon et al. 2007; 
Mannucci et al. 2008; Sand et al. 2011). 

The SN la number per formed stellar mass, A*'ia/M, can 
easily be related to the fraction, 77, of stars in some initial 
mass range [7711,7712], that eventually explode as SNe la: 

_ ^la Io°i m{dN/dm)dm 

M j2"^{dN/dm)dm ' ^ ' 

where dN/dm is the IMF. For the "diet Salpeter" IMF 
(which, again, gives results similar to the Kroupa IMF as- 
sumed by the VESPA SFH reconstruction), and an ini- 
tial mass range of 3-8 M©, often considered for the pri- 
mary stars of SN la progenitor systems, the ratio of the 
two integrals equals 33. Our present result of A^ia/M = 
0.00130 ± 0.00015 Mgi translates to 77 = 4.3 ± 0.5%. In- 
terestingly, Mannucci et al. (2006) estimated rj = 4.3% by 
modeling SN la rates per unit mass as a function of galaxy 
colors in a local sample. 

4.3 Analysis of cross-talk between bins 

In the Mil analysis of the LOSS SN sample and the VESPA 
SFHs of the monitored LOSS galaxies, a major system- 
atic error emerged, due to the limited 3" aperture size of 
the SDSS fibres. The LOSS galaxies are relatively nearby, 
z < 0.05, with half of them at 2 < 0.022, and are hence 



large in angle. The SDSS spectra in many cases are therefore 
dominated by old populations in the centres of the galaxies 
probed by the fibres. However, many of the LOSS SNe ex- 
ploded in the outer regions, where there is ongoing star for- 
mation outside the aperture of the SDSS spectroscopy. As 
a result, the DTD solution mistakenly associated these SNe 
with an old population, and hence a large delay. The result 
was an artificially enhanced DTD amplitude in the interme- 
diate and delayed bins, at the expense of the prompt bin, 
an effect referred to as DTD inter-bin "leak" or "cross-talk" 
in Mil. To mitigate this effect. Mil defined a culled sam- 
ple of galaxies that excluded Sab-Sbc Hubble types, which 
are most susceptible to this stellar-population misrepresen- 
tation by the fibre spectrum. The DTD derived from the 
culled galaxy sample, and the SNe that they hosted, had 
a much-reduced DTD leak, but at the price of a reduced 
sample of SNe, resulting in larger DTD Poisson errors. 

The DTD inter-bin leak, and its reduction in the culled 
sample, was diagnosed in Mil by deriving the DTD for 
core-collapse (CC) SNe. Core-collapse SNe explode within 
< 40 Myr of star formation. The CC SN DTD should have 
zero amplitude on timescales much longer than this, partic- 
ularly in the bins ai t > 420 Myr. Thus, any signal in those 
bins is due to a leak. Compared to LOSS, we expect the 
SDSS2 sample to have smaller or no inter-bin DTD leaks. 
The typical redshift of the SN hosts is z = 0.15, an order 
of magnitude larger than for LOSS, and all but three of the 
SDSS2 SN hosts are at z > 0.05. The spectra from the SDSS 
fibres should thus represent well the entire stellar population 
of each galaxy. To verify this, we have compiled a sample of 
all SNe in the SDSS2 website's list that have been spectro- 
scopically classified as CC SNe and are associated with a 
Stripe 82 galaxy with a VESPA SFH. There are 41 such CC 
SNe. This list is, of course, far from complete, particularly 
because the main SDSS2 science objective was the discovery 
of SNe la, and events suspected to be CC SNe had low pri- 
ority for spectroscopic foUowup (DIO). Nevertheless, even in 
a small and incomplete sample, DTD leak will still manifest 
itself in the form of a DTD signal at long delays. 

The leak effect can be quantified via the ratio between 
the intermediate and the prompt CC-SN DTD components, 
Vl'2/vl'i, which should be zero when there is no leak of signal 
from the first to the second DTD bin. In the full LOSS CC 
sample of Mil (see their table 1), *2/*i = 0.165 ± 0.045, 
i.e., there is clear presence of a leak. In the Mil sub- 



© 0000 RAS, MNRAS 000, 000-000 



SDSS2 Supernova Delay Time Distribution 11 



sample that was culled of Sab-Sbc Hubble types, this ra- 
tio was reduced to Vl'2/'I'i = 0.077 ± 0.050, which is con- 
sistent with zero. For the SDSS2 CC-SN sample, we find 
*2/*i = 0.069 ± 0.050, similar to the ratio in the culled 
LOSS sample, and again consistent with zero. 

As another test of cross-talk, we have recalculated the 
SN la DTD, limiting the galaxies and the SNe to be above 
some redshift threshold. Since the fibre-related DTD leak 
decreases with increasing galaxy redshift (as more of the 
galaxy's extent is included in the aperture), we expect a 
change in the DTD when changing thresholds, as was indeed 
seen in Mil. In the SDSS2 sample, however, we see no signif- 
icant changes in the SN la DTD when limiting the sample 
to higher redshifts. We conclude that the SDSS2 sample's 
DTD is largely free of the inter-bin cross talk due to the 
SDSS fibre aperture. 

4.4 Stretch-dependent DTD 

Most DTD recovery work to date has considered only a uni- 
variate, universal, DTD. A correlation has been known for 
some time (e.g. Hamuy et al. 2000; Neill et al. 2009; Hicken 
et al. 2009) between the age of a SN la host and the SN lumi- 
nosity (or, equivalently, its light curve's "stretch parameter", 
s, or the mass of radioactive Ni synthesized). Old galax- 
ies tend to host low-luminosity SNe la, while star-forming 
galaxies tend to have luminous SNe la. This is an observed 
link between the progenitor population and the energy of 
the explosion, and thus may provide a critical clue to the 
progenitor question. However, it is only through the DTD 
recovery approach that one can distill this link from the 
data, decoupled from the SFHs of the individual galaxies. 
BIO used the stretches of their SN la sample's light curves 
to divide their sample into high-stretch and low-stretch sub- 
samples, and derived the DTD for each subsample. This 
was the first derivation of a bivariate '^{t,a) distribution, 
one with three delay-time bins and two stretch bins. Here, 
"DTD" is no longer an appropriate name, as this is now the 
bivariate distribution (or "response function" ) of delay times 
and stretches. The bivariate response contains information 
that is additional to the distribution's univariate projection, 
the DTD, as it gives not only the age of the progenitor sys- 
tems but also the run of explosion energies for each progen- 
itor age. BIO indeed found that luminous, high-stretch, SNe 
la tend to have most of their DTD power at short delays, 
while low-stretch, underluminous, SNe la have a DTD that 
peaks in the longest-delay bin. 

To test this result with our newly defined SDSS2 sample 
and analysis, we have created low-stretch and high-stretch 
SN la subsamples, adopting, like BIO, a border between the 
two populations of s = 0.92. Among the 61 of our SNe that 
are from the DIO sample with spectroscopic confirmation, 
50 are in the BIO sample, and we adopt the BIO classifica- 
tions for these SNe as high- or low-stretch. We exclude the 
remaining 11 DIO SNe from the present analysis. For the 
62 photometrically classified SNe la that are from Sll, we 
convert the AM15 values, tabulated by Sll, to stretch, s, 
using (Guy et al. 2005, their footnote 7), 

s = -0.59 AA/i5 + 1.55. (12) 

The 9 additional SNe from our full sample, that were pho- 
tometrically classified by DIO, do not have stretch, or anal- 



ogous, parameters tabulated, so we also exclude them from 
the present analysis. With this selection, we obtain 45 low- 
stretch SNe la and 67 high-stretch SNe. 

Table [2] lists the DTDs we obtain separately for the 
low-stretch and high-stretch subsamples. The high-stretch 
DTD has a 3.5cr signal in the prompt time bin, and ~ 2.5cr 
detections in the later time bins. The best-fitting power law 
to this DTD has an index -1.12±0.13. The low-stretch DTD 
has only 1.5(t detections in the prompt and intermediate 
time bins. Those DTD components are thus consistent with 
zero, but higher ^'i and ^'2 values, giving similar component 
ratios as in the high-stretch DTD, also cannot be ruled out. 
The delayed low-stretch component, '^■i, does have a > So- 
signal. A power-law index in this case is — 0.9^015. 

Qualitatively, these results are similar to those of BIO, 
in the sense that the high-stretch SNe have a DTD with 
a strong prompt signal, while the low-stretch SNe have a 
strong signal in the delayed bin. However, the present analy- 
sis does not reproduce the nearly one-to-one correspondence 
reported by BIO, where essentially all high-stretch SNe de- 
rive from a young population, and all low-stretch SNe from 
an old one (see also further discussion of the BIO results in 
Section [4.51 below). Rather, in the present analysis, there 
may be a trend in the above sense, but high-stretch SNe 
do have significant non-zero amplitudes in their interme- 
diate and delayed DTD bins. And, within the current un- 
certainties, low-stretch SNe may have non-zero prompt and 
intermediate-delay DTD components. 

We suspect that the difference between the present re- 
sults and those of BIO may be due to selection effects that 
arise from the use, by BIO, of a sample that includes only 
spectroscopically confirmed SNe la. The SDSS2 survey was 
strongly oriented toward dicovering SNe la for cosmographic 
purposes. The selection of candidates for spectroscopic fol- 
lowup, and the success of those followup observations, were 
both likely tied to the luminosity of the SN (and hence to 
its stretch) and to the type of host galaxy. For example, a 
candidate in an early-type galaxy would have been likely 
considered a promising SN la, even if its luminosity was not 
high. Its spectroscopy may have been more likely to succ- 
ceed in the presence of the lower host starlight background, 
compared to a similar low-luminosity candidate in a star- 
forming galaxy. Such biases could artificially eliminate, in 
the spectroscopically confirmed sample, low-luminosity SNe 
la in star-forming galaxies, and this would reduce or elimi- 
nate the ^1 component in the low-stretch sample. The more- 
complete, spectroscopic plus photometric, sample analyzed 
here is more immune to such a possible selection effect, and 
may thus be reflecting better the true DTDs. Recently, Per- 
rett et al. (2012) have also failed to find strong differences in 
volumetric rate evolution between low-s and high-s SNe la. 
Larger, well-defined, samples will be able to further probe 
the bivariate distribution of SN la delay and stretch. 

There is some concern that the DTD that we have recov- 
ered may be biased, due to the known correlation between 
host galaxy age, or star-formation rate, on the one hand, 
and SN stretch or luminosity, on the other. The DIO detec- 
tion efficiency function that we have used does not account 
for this correlation. Thus, for example, early-type galaxies, 
which tend to host lower-luminosity SNe, would have a true 
detection efficiency lower than that of DIO. Our overestimate 
of the efficiency could then lead to a perceived paucity of SNe 
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la in old populations, and hence to an underestimate of the 
true DTD component. To gauge the possible magnitude 
of such an effect, we have repeated the DTD derivation, but 
modifying the DIO detection efficiency function as follows. 
Sullivan et al. (2006) have shown (their fig. f 1) that galax- 
ies with specific star-formation rates sSFR > fQ~®'^yr^^ 
host SNe with a mean stretch of s ~ i.03, while the most 
passive galaxies, with sSFR < 10^^'^yr^^ have SNe with 
mean s « 0.87. This change in s corresponds to a change 
in SN peak absolute magnitude by about 0.25 mag. Such 
a change in brightness, in turn, corresponds to a change in 
redshift of Az/z = ±6%. For galaxies that, based on their 
VESPA SFHs, have current specific star-formation rates 
sSFR > fO~^'^yr~^, we therefore modify the DIO efficiency 
function (described in Section [2. such that it declines lin- 
early from z = 0.184 (rather than z = 0.175), reaching zero 
at 2 = 0.42 (rather than z = 0.4). For the passive galax- 
ies with sSFR < lO^^^yr^^, we assume a detection func- 
tion that declines linearly from z — 0.166, reaching zero al- 
ready ai z — 0.38. For the intermediate, weakly star-forming 
galaxies, we use the DIO detection efficiency, as before. We 
find that the recovered DTD is only slightly affected by this 
modification of the detection efficiency, with ^'i decreasing 
by 10% and increasing by 10%. The change leads to a 
slightly less steep power- law fit to the DTD, with an index 
oi 13 = -1.07 ±0.07. 

4.5 Comparison of the BIO and Mil DTD 
recovery methods 

As elaborated above, the SN la sample and the details of 
the DTD calculation based on SDSS2 SNe differ in a num- 
ber of ways between the present work and that of BIO. Nev- 
ertheless, it is interesting to compare the performance and 
results of the two different algorithmic approaches - on the 
one hand, that of Mil and of this work, i.e., comparing ob- 
served and predicted (given a DTD model) SN numbers in 
each galaxy; and on the other hand that of BIO, compar- 
ing the observed mean host spectrum to a mock mean host 
spectrum (for a given DTD model). We have therefore cal- 
culated a DTD reconstruction with the Mil approach, but 
reproducing as closely as possible the SN la sample and the 
assumptions of BIO. 

To construct a SN la sample analogous to the one of 
BIO, we have matched coordinates and redshifts between 
the 101 SNe la in table 2 of BIO and the current Stripe 82 
VESPA sample of galaxies (without any RA limits). We have 
matched 89 of the BIO SNe to host galaxies (we note that 37 
of the SNe in BIO have mis- typed declination signs - positive 
instead of negative). This smaller number of SNe, 89 com- 
pared to 101, is expected, given our smaller galaxy sample 
(66,000) compared to that of BIO (77,000). Among these 89 
SNe la, BIO listed 32 as low-stretch and 57 as high-stretch. 
We have then run these low-stretch and high-stretch samples 
through the Mil DTD recovery algorithm, but using the 
same efficiency function assumed by BIO, with the same pa- 
rameters that they assumed for high-stretch and low-stretch 
SNe. In this calculation, we do not correct for cosmological 
time dilation, and do not scale down the VESPA masses by 
0.55 (see §2, above). We find for the three of components 
of the DTD, *i, *2, and *3, respectively, in units of 10"" 



Low-stretch: (0 ± 28; ± 5; 1.3 ±0.4); 
High-stretch: (65 ± 15; 8 ± 3; 0.1 ±0.2). 
Table 3 and figure 7 of BIO give, for several of the low- 
stretch and high-stretch ^ components (eh,2, and ei,2, in 
their nomenclature) their results and uncertainties. To these 
uncertainties we have added la Poisson errors, not included 
in BIO. The BIO *i results are. 
Low-stretch: (0 ± 20; ± 5; 1.6 ±0.3); 
High-stretch: (75 ± 19; 8 ± 5; 0.3 ±0.2). 
We see that the best-fit results and the uncertainties found 
by the two algorithms are quite similar. The slightly higher 
values found by BIO for the significantly non-zero compo- 
nents may be due to the fact that BIO draw the DTD values 
that they test from an assumed prior distribution, while in 
Mil and here we do an adaptive grid search for the best 
fitting DTD in the 3D space of *i, *2, and *3. 



5 DISCUSSION AND SUMMARY 

We have taken the DTD recovery method presented in Mil, 
where it was applied to the LOSS SN sample, and applied 
it here to a subsample of SDSS2 SNe la that were hosted by 
Stripe-82 SDSS galaxies with VESPA SFH reconstructions. 
Among all of the DTD recontructions to date that are based 
on detailed SFHs of individual galaxies (BIO; Mil; Maoz & 
Badenes 2010; this work), we have obtained the DTD that 
is the most precise (in terms of Poisson errors) and the most 
accurate (in terms of control of systematic errors). The full 
Mil sample was systematically affected by the limited SDSS 
fibre aperture, which led to cross talk between DTD bins. In 
the reduced Mil sample, where this problem was mitigated, 
only 49 SNe la remained, leading to large Poisson errors. In 
the analysis by Maoz & Badenes (2010) of 77 Magellanic 
Cloud SN remnants, the bulk of the DTD signal was at 
t < 35 Myr, corresponding to core-collapse SNe. This meant 
that only of order 10 SNe la were effectively driving the SN 
la DTD, which was therefore very noisy. Another possible 
systematic error in Maoz & Badenes (2010) was contamina- 
tion of the SNa la DTD by CC SNe, with progenitors near 
the poorly known progenitor age border between CC-SNe 
and WD formation. 

The SNe from SDSS2, which are at larger distances and 
are well classified, can avoid the above systematics, and at 
the same time allow for relatively large numbers of SNe la. 
The BIO analysis of SDSS2 SNe, however, used a limited 
sample, available at the time, consisting only of spectrally 
confirmed SNe. Furthermore, BIO estimated a detection effi- 
ciency function that did not take into account the details of 
the real observing patterns, the detection process, classifica- 
tion, and sample selection, details that were later quantified 
by DIO, based on real-time detection simulations. Finally, 
the BIO analysis omitted necessary corrections to the SN 
visibility times and to the stellar masses of the galaxies. 
All of the above are potential sources of systematic error, 
affecting each DTD component by up to a factor ~ 2. In 
our renewed definition and analysis of an SDSS2 SN sam- 
ple, we have addressed these issues. In this analysis, we have 
detected for the first time, at a highly significant level, an 
intermediate-delay DTD component, at 0.42 < t < 2.4 Gyr. 
BIO and Mil observed such a component only at the ~ 2a 
level. In addition, the completeness of the sample and the 
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consideration of absolute stellar mass, visibility time, and 
efficiency, permit a robust derivation of absolute DTD lev- 
els. 

It is interesting to consider the cause behind the much 
more significant detection of an intermediate-delay DTD 
component, ^'2, compared to the BIO analysis of a differ- 
ently defined SDSS2 SN la sample. Given that the Mf f and 
BIO algorithms give similar results and uncertainties (Sec- 
tion [43] above), the difference in the ^2 significance could 
be due to the differences in the SN samples, the differences 
in the adopted efficiency functions, or both. To investigate 
this, we have recovered best-fit DTDs with the present SN la 
sample (Table f), but using the BIO efficiency function, and 
omitting other corrections, as in Section [4.51 Alternatively, 
we have used the BfO SN sample as described in Section Pl.Sl 
but using the empirical DIG efficiency function, 1 + z correc- 
tion, and so on, as in our main calculation. We find that nei- 
ther choice, alone, raises much the ^'2 detection significance. 
Apparently, the improvement results from the combination 
of a carefully selected, relatively complete, sample of SNe, 
and a realistic efficiency function, one that embodies the se- 
lection effects of that same SN sample. The two, together, 
sharpen the correlation between the presence of SNe la and 
an intermediate-age stellar population in some galaxies, and 
thus reveal clearly the ^'2 component. 

Our findings point to a continuous distribution of SN 
la delay times. Although recent observational developments 
have led to discussions of "prompt and delayed SNe la", 
"two SN la channels/populations", "a bimodal DTD", and 
so on, theoretical binary population synthesis models have 
almost always predicted a broad and continuous DTD, even 
for a single physical channel (e.g. DD). Thus, the continuous 
observational DTDs that are emerging (Totani et al. 2008; 
Maoz et al. 2010; this work) are not unexpected. Naturally, 
DTDs recovered in the future with finer time resolution may 
still reveal a real bimodality. The current coarse time reso- 
lution of the DTD is dictated by the resolution of the galaxy 
SFHs. Future improved data will also permit recovery of a 
better-resolved, and higher signal-to-noise, bivariate distri- 
bution of delay and stretch, ^'(i, s), that could sharpen the 
current hints of a relation between the progenitor age and 
the explosion energy. 

As in most of the recent studies, the DTD we find ap- 
pears to be consistent with the ~ form generally ex- 
pected in the DD scenario. This result joins others that 
are strongly suggestive of DD progenitors for at least some 
SNe la. Most notable among these are the absence of any 
signs of the putative donor star, a few years before, and dur- 
ing the explosion of SN2011fe (Li et al. 2011c, Nugent et al. 
2011, Horesh et al. 2012), and 400 years after the explosion 
that left the SN remnant 0509—67.5 (Schaefer & Pagnotta 
2012). It is noteworthy that the main objection to the DD 
model has been theoretical, arguing that a DD merger can- 
not produce an event with the characteristics of SNe la. 
The fact that there are, by now, a couple of cases where 
we may have witnessed DD mergers producing normal SNe 
la suggests that the theoretical objection could be somehow 
fiawed (unless the absence of a companion, in these cases, is 
the result of one of the "SN on hold" scenarios mentioned in 
Section 1, assuming they work). If indeed the DD scenario 
is physically possible after all, perhaps most or even all SNe 
la could result from DD mergers. Examining the evidence 



from the opposite direction, Maoz et al. (2012) and Badenes 
& Maoz (2012) have recently estimated the specific merger 
rate of Galactic WD binaries. They found that, if one consid- 
ers all WD mergers, most of which have sub-Chandrasekhar 
merged masses, but only a few tenths of a solar mass below 
this limit, their merger rate is remarkably similar to the spe- 
cific SN la rate of Milky- Way-like galaxies. Naturally, this 
coincidence is meaningful only if sub-Chandrasekhar CO- 
CO WD mergers can indeed produce a SN la, which is still 
an open question (e.g. van Kerkwijk et al. 2010; Guillochon 
et al. 2010; Pakmor et al. 2011). 

Another interesting result of our analysis is the absolute 
DTD normalization that we find - lower than most of the re- 
cent determinations based on SN la rates in nearby galaxies 
and in galaxy clusters, but consistent with the levels found 
based on volumetric SN la rates. This result relaxes some 
of the tension between the previous, high, observational SN 
la production efficiencies, and the lower predictions from bi- 
nary population synthesis models (see, e.g. Maoz 2008; Men- 
nekens et al. 2010). On the other hand, our result re-floats 
the possibility that SNe la, for some reason, are produced 
more efficiently in galaxy cluster environments than in fleld 
environments. If this is true, it could be the result of an 
IMF that depends on environment. Evidence has accumu- 
lated recently that the most massive elliptical galaxies have 
"bottom-heavy" IMFs, with low-mass components as steep, 
or even steeper than, the Salpeter IMF (e.g. Treu et al. 2010; 
Cappellari et al. 2012; Conroy & van Dokkum 2012). Clus- 
ters host the most massive ellipticals. Thompson (2011) has 
proposed that a large fraction of SN la progenitor systems 
are triples, in which a low-mass tertiary induces ellipticity 
in the inner WD binary via the Kozai (1962) mechanism, 
signiflcantly increasing the time-integrated number of DD 
systems that merge (and presumably produce SNe la). A 
speculative possibility that we raise is that, in cluster galax- 
ies, the large availability of low-mass stars translates into 
more tertiaries, and hence a larger SN la production effi- 
ciency. 

To summarize, we have constructed a new SDSS2 SN la 
sample that is more complete and has a better-characterized 
detection efficiency than what was available in BIO. We have 
analyzed it with the methodology of Mil, but avoiding the 
small sample size and the systematic errors in SFH that 
were unavoidable in the LOSS sample of Mil, while making 
several necessary corrections that were omitted in BIO. In 
the end, we find a SN la DTD that is more accurate and 
more precise than the ones found by either BIO and Mil. 
For the first time, a DTD based on the three time bins used 
here shows a highly significant signal in all three bins. The 
DTD is consistent in form with the found in the stud- 
ies reviewed above (including BIO and Mil), but differs in 
absolute normalization from some of those studies. 

The SN la progenitor problem is far from solved. How- 
ever, our results demonstrate how our view of the SN la 
DTD can progressively improve by means of analysis of sam- 
ples of SNe that are progressively larger, purer, more com- 
plete, and well-characterized, and with samples of monitored 
galaxies with SFHs that are better determined, systemati- 
cally, temporally, and spatially. These improved views will 
continue to play a role in resolving the progenitor issue. 
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